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THE SOLAR REFLECTANCE OF A SNOW FIELD 

B. J. Choudhury 
Computer Sciences Corporation 

A. T. C. Chang 

NASA, Goddard Space Flight Center 
ABSTRACT 

The radiative transfer equation has been solved using a modified 
Schuster-Schwartzschild approximation to obtain an expression 
for the solar reflectance of a snow field. The parameters in the 
reflectance formula arc the single scattering albedo and the fraction 
of energy scattered in the backward direction. The single scattering 
albedo is calculated from the crystal size using a geometrical optics 
formula and the fraction of energy scattered in the backward direc- 
tion is calculated from the Mie scattering theory. Numerical results 
for reflectance are obtained for visible and near infrared radiation 
for different snow conditions. Good agreement has been found 
with the whole spectral range. The calculation also shows the 
observed effect of aging on the snow reflectance. 
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THE SOLAR REFLECTANCE OF A SNOW FIELD 


INTRODUCTION 

Inica'st in developing techniques for remote measurement of snow parameters (e.g. 
temperature profile, density grain size, water equivalent and free water content), have 
significantly increased in recent years. One of the more promising techniques to monitor 
these parameters is that of microwave radiometry. A recent study by Chang et al.. (1976) 
has used a theoretical scattering model to relate the observed microwave brightness tem- 
perature to the physical temperature and snow grain size for snow on glacier ice. By 
utilizing this remotely determined snow grain size and density information, the solar 
spectral reflectance from a snow field can be Inferred. The knowledge of spectral reflec- 
tance is important because it determines the amount of solar energy that will be retlected 
by a snow field. The present study is an attempt to calculate the solar reflectance of a 
snow cover in terms of parameters which depend upon the grain size and density. 

For incident solar radiation, the energy reflected at the surface constitutes only a small 
fraction (about .** percent) of the total energy that is reflected by a snow field. For 
visible and near infrared radiation, the scattering by individual snow crystals is sucli that 
major portion of the radiation gets scattered in the forward direction. The small traction 
of energy that is scattered in the backward direction, however, gets enhanced considerably 
by the cumulative effect of many scatterers. Thus the major portion of the incident solar 
radiation which gets reflected by a snow cover is due to the volumetric effect of many 
scatterers. 
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The appropriate equation which needs to be solved for calculating the cumulative effect 
of many scatterers is the radiative transfer equation. An analytical solution of the 
radiative transfer equation when the individual scatterings are in the near forward direc- 
tion is rather difficult, but a solution can be obtained by numerical methods (Choudhury, 
1977). In this paper we will discuss an analytical solution with modified Schuster- 
Schwartzschild approximation (Sobolev, 1963; Sagan and Pollack, 1967). The expres- 
sion for spectral reflectance obtained from this analytical solution of the radiative 
transfer equation can be evaluated in terms of snow parameters. The theory presented 
here is intended to show the influence of snow parameters on the solar reflectance of 
snow. 

In the following the general expressions and equations needed to study the spectral reflec- 
tance of a snow cover are given. It also includes the background information regarding 
the snow parameters and the solar intensity distribution which are expected to influence 
- the spectral reflectance. Tlie analytic solution of the radiative transfer equation will be 
discussed and an expression for the spectral reflectance obtained. In addition comparisons 
are made for the obtained numerical results and the observations. 

ulNeral formulation 

The solar energy incident on the snow surface can be divided into two components: 
the collimated intensity. 1^, falling at an angle 0 q with respect to the normal and the 
nearly isotropic diffuse intensity 1^ which arises due to the atmospheric scattering of 
the radiation. The tlux of the incident radiation is given by 



P| " *c *d » 


( 1 ) 


where Mq * cos 6 q . 

Since most of the incident radiation is contained in the collimated component, the flux 
of the radiation reflected by the snow-air interface is approximately given by 

Fr = r(0o)dc#^o *d (2) 

where r ( 6 ^ is the Fresnel reflectivity of the surface. 


The refracted portion of the incident radiation will progress through the medium (snow) 
getting weaker due to scattering and absorption. It is convenient to separate the refracted 
portion of the collimated radiation into two components; the unattenuated radiation and 
the radiation which arises due to scattering. This separation is done because the unatten- 
uated radiation remains angularly concentrated along the direction of refraction but the 
angular distribution of the scattered radiation depends upon the scattering phase function. 
If multiple reflection at the boundaries of the medium is negligible then the unattenuated 
intensity at depth x is given by: 


lu (’■) = (1 - r(flo)) n^ Ig 5 (0 - 0q) 6 (// - Mq) exp (-r//iy) 
where n is the index of refraction, /i|j is the angle of refraction: 

^0= » 


(3) 


. izA 

«2 


'/i 


and r is the optical depth 


■f 

Jo 


Te(x') dx' 


(4) 
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when* Y^.(x) is the extinction coefficient at depth x. The delta functions in (3> repre- 
sent the constraint due to the collimated nature of the radiation. 


The t1ux of the unatteiuiated radiation is given by 


I ’u ' ( • - r Ic ^o 


(5) 


riic intensity which arises due to scattering of the incident radiation is to be obtained by 
striving the radiative transfer equation (('handrasekhar. I‘*o0): 


dl (r^i) 

H ^ = - l^lT.^i) 

dr 


^ j P tp.Ai') Ij (r./u'> d*i' 




exp 


(M 


where n is the angle of radiation, u? is the single scattering albedo and p(^./i’l is the 
scattering phase function. 


The boumlary conditions for the radiative transfer equation are as follows: 

fl) If the snow-air interface is assumed to be a geoinetncally smooth surface then 
the appropriate bouiulary condition is 

li (r = 0. ^) = r (^H_ (r = 0. + ^l-r(^>)n* l^j (7) 

where 4 is tlie radiation going downward into the snow and l_ is the radiation going 
upward towaiils the atmospheie I'lie first term on the right hand side represents the 
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intensity reflected by the snow-air interface and the second term accounts for the re- 
fracted component of the diffuse intensity. 

(2) The second boundary condition is associated with the snow-soil interface. 
The general boundary condition will be of the form: 

L (Tq, fi) = in) 4 (t^,. n) (8) 

where r, in) is the reflectivity of the snow-soil interface and Tq is the total optical 
thickness of snow. 


Once the solution of equation (6) is known subject to the above boundary conditions, 
one can calculate the flux of the radiation emerging from snow as: 


-i: 


L(o./i) (i - Tin)) 


ndn 


i')) 


The albedo of snow is to be calculated using equations (1), (2) and (9) as: 


A = 


F: 


( 10 ) 


The volumetric effect of scattering which is mainly responsible for large values of albedo 
is associated with the emergent flux appearing in equation (10) and is defined through 
the equation (9). The calculation of the emergent flux requires the solution of the 
radiative transfer equation (6). 


AP PROXIMATE EXPRESSION FOR SNOW REFLECTANC E 

To obtain a closed form analytic solution of the radiative transfer equation (6), the 
following assumptions are made. 
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1. The medium (snow) is u homogeneous plane parallel layer, consisting of 
statistically equal size Ice spheres. 

2. The Schuster-Schwartzchild or two-stream approximation for the intensity is 

valid. 

The angular asymmetry in the scattering process can be described by introduc- 
ing a parameter which weights the scattered intensities. 

It should be noted that with the first approximation we aa* neglecting the presence of 
all internal inhomogenities like enlarged ice glands and ice lenses which sometimes form 
during the process of luetamorphism. The other two approximations are assticiated 
directly with the method of solution. It has been shown that for a optically thick 
medium, the solution of the radiative transfer equation under the above set of approxi- 
mations is reasonably accurate (Sagan and Pollack 1^)07). 

Hie Schuster-Schwartzschild or twivstream approximation is based upon partitioning the 
scattered radiation into mean intensities along the forward and the backward hemispheres. 
In single scattering, if 0 is the fraction of energy scattered in the backward hemisphere 
with respect to the direction of incidence, then we can replace the radiative transfer 
ciiuation ((') In the following equations 

d4 

= - Ia + UJ( 1 - d) 1+ W>d 1 
dr 

10 ( 1 - d) I' exp (- r) (11) 


h 


and 


- - I ♦ w(I L ♦ CO/JI+ + 

dr 

W0I' exp t-T) ( 12) 

where 

r = I I -r(o)| (I,. + Ij) 

and I_ arc respectively the mean intensities along the forward and backward direc- 
tions with respect to the direction of incident radiation, and r(o) is the Fresnel retlec 
tivity for normal incidence. It should be noted that these equations represent the 
radiative transfer in a one dimensional medium (Sobolev. I^h>.t). 

The radiative transfer equation (1 1) and (12) contains two parameters. nanieU. (I) cc. 
the single scattering albedo and (2) (i. the fraction of energy scattered in the backward 
direction. We will discuss the metluHl of determining these parameters before striving 
this radiative transfer equation. 

For ice particles in the inilliinetcr range and the radiation in the visible or near infrared 
legion, the si/e of the particle is significantly larger than the wavelength. I'lulei this 
condition. Mie scattering theory shows the extinction cross-sect it'ii undergoes resonance 
oscillations for small changes in the wavelength and size of particle. This restMiance 
elTect, however, can be blurred by assuming tlut there is a statistical variation in the 
particle size about its mean rad«u». With this mean size value, one can use the geomet- 
rical optics approximation to calculate the single scattering albedo using the equation 
(Irvim ,ind Pollack. I*>b8. Sagan and Pollack, 
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tJ = 1/2 + 1, 2 • exp (- 2kj^ r) 


(I.^) 


where r is the radius of the particle and the absorption per o.iit length of icc at 
wavelength. X. In our calculation, we have used kj^ values tabulated by Inine and 
Pollack (l^oNI. I'ollo ving above discussion, instead of using equation (13). it should 
be noted that one can also us«* Mic scattering tiuorx to calculate the single scattering 
albedo. 

(or.sistent with equation (l.t) the extinction a>efficient required in the calculation of 
'•[ tical depth (equation 4) is given b> 


>e 



( 14i 


where pj. and pj are res|H*i'tiveIy the density of snow and the density of ice. 


The fraction of energy scat:ercd in the backward direction, d. in a single scattering has 
been studied b; oiany inve (e.g. Barlky. I'^t'S. S.igan and Pollack. l%7i. This 

parameter can be related to the scattering phase function p(cosO) as. 


( cos () ) • cos 0 7~ ( 1 .s ) 

47T 

where 0 is the scariering angle and dS2 is an element of solid angle in the ci'-ordm.ite 
system wtiere ft is the pt'la. angle 1 he norniali/ation condition for the phasi* function 
is 



Jp (ct>s(M 


dfZ 

4ir 


1 . 


S 


When the phase function is expanded in terms of Legendre polynomials, the parameter 
P can then be expressed in terms of the coefficient of the first order Legendre poly- 
nomial. 


For a given particle size one can use Mie scattering theory to calculate the phase func- 
tion and then use (IS) to calculate the value of p. Calculations performed by Irvine 
and Pollack ( 1968) show that when the size of the particle is significantly larger than 
the wavelength of the radiation, the value of P is about 0.07S, i.e., in single scattering 
only 7 percent of the incident energy is scattered in the backward direction. In our 
calculation, we have used this value of P (0.075) because for visible and near infrared 
radiations, the scattering by individual ice crystals is indeed in the near-forward direction. 


For a thick snow layer, the general solution of equations (11) and (12) can be obtained 
by standard methods (see Sobolev, 1963 pp. 33) as 

4 (t) = 1/2 (1 + *-^) C e'^’^ - r e*^ (16) 

and 

I. (T) = 1/2 (1 - -^^) Ce**^’’ (17) 

where C is the constant of integration and 

k^ =(1 -w) (l -w(l -2/3)) 


To determine the constant of integration we now impose the boundary' condition (7) to 
obtain 


c = 


((’• 


2J^ 

1 - (jJ 1 - CJ 

+ — r- )-r(o) (1 --T— 


)) 


(18) 
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With the knowledge of ccunplete solution, the reflectance for a onc-ilimensional system 
can bo calculated as: 


A ■ r(o) + l_ tr » o) (l - rtol) ^ 


The explicit expression for reflectance is: 


(d - r(o)) ^ tt 
A » r(o) + ' — . J 


where 


- rto)a 


a 


k - I + oj 

k + 1 - to 


(Note that for u? = 1 one finds a = 1) 


tl'J) 


(201 


I hc 1‘resncl reflectivity for normal incidence. r(o). can be calculated from the reflective 
iiulex of ice. n. using the formula; 


r(o) = 


n - I 
n I 


cn 


.Since the ^|vctral dc|HMidence of n in the visible and the near infrared region is extremely 
weak (Irvine ami Pollack. I'^oS). we have used a constant I..U) for n in calculating the 
rtol. 


The expression for reflectance is in close n'seinblance with that calculated by Dunkle and 
bevans tl'J.'vM without explicitly making reference to one-dimensional system. The ilefini- 
lion of the parameters appearing in evpiation (20) are. however, quite different. We thus 
luvte that although theie Is a formal resemblance between the expression for refleclana' 
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obtained in this paper and that obtained by Dunkle and Revans (1Q5o). the method of 
obtaining this expression and the nrethod of calculating the parameters of this expression 
are quite different in two cases. 

The expression for reflectance obtained in this pa|ier is based upon simplifying approxi- 
mations. We have discussed these approximations and attempted to provide their justifi- 
cation. Tire ultimate justification of the model is to be found by aunixiring with the 
experimental obserx-ations. In the next section, we will provide the numerical results 
and compare them with the observations. 

? ANjl DISC USSIONS 

By examining equation 20. the single scattering albedo u is tlie only varx ing parameter 
in calculating the snow reflectance for a given wavelength. Since co is directly n.'lated 
to the snow crystal radius, it is pc*ssible to char.icteri/e the snow condition I'y s|H*cifying 
the radius of the ice crystals. 

file ice crystals of fresh fallen snow usually are of complex form aiul contain sharp 
corners. I'he shape of these crystals changes with time de|KMuling upon the vapor cim- 
tent and the prevalent temperature. The process of equi-temperatiire metamorphisni 
leads to the production of fairly uniform and well rounded grains. At the initial stage 
of this metamorphism the mean radius of the icc crystals is about 0.2 mm. Tlie radius 
then c^ontinues to incrc.ase as the process of metamorphism advances. At a fairly ad- 
vanced stage of this metamorphism, the radius increases to about 1.0 mm. I'he effect 
of other metamorphism such as the temperature-gradient metamorphism or the 
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mcIt-fivc /0 melamorphism. is itcnorally to pn>duv*c non-sphcrical and non-uniform ice 
crystals also the crystals arc laixcr than CMui-tcm|\'raturc mctanior|>hism. Thus, althoutih 
it IS not unitpic, one should he able to use the radius of ice crystal to characterize the 
si;g!c of inctamorphisin. 

l-iisure I illustrates the changes in the calculated snow retlectance due to the difference 
in the crystal radius. The chosen radius are the typical snow crystal sizes for different 
stajjcs of the enui-tem|H'ratuie inetamorphism. The overall shape of the curve diX’s not 
seem to deiXMul upon the railius of the ice crystals but the relative magnitudes of various 
parts of the curve var> with the radius. I he spectral ivtleclance in the red and near- 
infran'd regions shows nia\iinuin sensitivity to the cr>-stal sizes. 

In I'igure 2 we show the calculated snow reflectance compared with the exjH'riiiKMital 
values for a nearh fresh snow (O’Hrien and Munis. The a'sults of the calculation 

are in gv'i'd .igreeinent with the observations. .Ml prominent sjx'ctral structures appearing 
in the observatum are well duplicated in the calculations. The ituantitative agreement is 
also good but should be treated with caution because the e\|X'riniental results are relative 
to a stand.ird (white Mariuni Siilfatet reflector. It is however n*assuring that the radius 
values used in the calculation are in the range expected for a fresh fallen snow. 

In l igure .' we compare the calculated and the observed reflectance of a naturally .igeil 
two days old snow (OMtiien aiul Munis. I'J'*.^). There is a goixl qualitative agreement 
between them The radius for the snow crystals which give agreement with the obsigna- 
tion are l.iigei than those for the fivsh fallen snow (I'igure 'T The natural aging prv'cess 


of snow, during the observations, was such that the ambient air temperature was hover- 
ing above and below freezing point. It is unlikely that the snow was undergoing the 
equi-temperature metamorphism. As a result, the ice crystals were probably non-uniform 
within snow. Furthermore, if the effect of the solar irradiance is interpreted as simply 
to advance the process of metamorphism of the snow layer in contact with the air, 
then the radius of the ice crystals in the top layer will be larger than those deep within. 
Since the absorption coefficient of ice increases with the wavelength, the effective thick- 
ness of the snow layer that contributes to the reflectance decreases with the wavelength. 
Therefore, the reflectance for longer wavelengths should have more influences from the 
top layer of snow and the reflectance for shorter wavelengths should correspond more to 
the deeper snow layer. This is a probable cause for the reflectance of the smaller crystals 
giving better agreement at shorter wavelengths while the reflectance of the larger crystals 
giving better agreement at longer wavelengths. 

Based on calculated results, two potential applications in remote snow reflectance 
monitoring are possible: (1) the use of multispectral bands to infer the radiation heat 
flux over a snow covered area; and (2) the use of visible and near-infrared channels to 
detect melting snow. 

The large contrast in snow and soil reflectance makes it possible to monitor quantitatively 
the amount of solar radiation absorbed by the earth’s surface. The degree of heat insula- 
tion by snow cover depends on the snow reflectance which is shown to be directly related 
to the snow crystal size in this paper. Therefore a spacebome multispectral scanning 
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instnmK'nt can provuic a method to monitor the heat tlux over a anow cover area 
under the National weather and ciimate program. 

The sea>iul potential application, that of detecting melting miow. is based on the calcula- 
tion that the snow reflectance decreases rapidly with the increase in snow crystal si/e. 
I'he existence of melt water on the snow surface will tend to increa.se the crystal si/e, 
also in the near infrari*d a'gion the absorption per unit length of water is higher than 
that of ice. As a result, the reflectance of melting snow is expected to be lower than 
the dry snow in the near infran'd region. 

I'he theory presented here, does show the qualitative and quantitative features of the 
ob.si'ncd reflectance. However the one dimeirsional approach can not distinguish the 
eflects of the direct and the diffuse solar radiations. \ more complexed analysis is 
needed to as.sess the elTect of directional irradiation. 
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FIGURE CAPTIONS 


Figure I. Illustrutions of the Effect of Different Snow Crystals on Snow Reflectance. 
Figure 2. Comparison of Calculated and Observed Reflectance of a Nearly Fresh Snow. 
Figure 3. Comparison of Calculated and Observed Reflectance of a Naturally Aged Snow. 
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